Variational approach to the modulational instability 
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We study the modulational stability of the nonlinear Schrodinger equation (NLS) using a time- 
dependent variational approach. Within this framework, we derive ordinary differential equations 
| (ODEs) for the time evolution of the amplitude and phase of modulational perturbations. Analyzing 

the ensuing ODEs, we re-derive the classical modulational instability criterion. The case (relevant 
to applications in optics and Bose-Einstein condensation) where the coefficients of the equation are 
, time-dependent, is also examined. 



Modulational instability (MI) is a general feature of discrete as well as continuum nonlinear wave equations. This 
^vj instability shows that in such settings, a specific range of wavenumbers of plane wave profiles of the form u(x, t) ~ 
exp(i(kx — u>t)) becomes unstable to modulations. The latter effect leads to an exponential growth of the unstable 
modes and eventually to derealization (upon excitation of such wavenumbers) in momentum space. That is, in turn, 
equivalent to localization in position space, and hence the formation of localized, coherent solitary wave structures 
[!]• 

The realizations of this instability span a diverse set of disciplines ranging from fluid dynamics [2] (where it is 
usually referred to as the Benjamin- Feir instability) and nonlinear optics [3] to plasma physics [4]. One of the earliest 
contexts in which its significance was appreciated was the linear stability analysis of deep water waves. It was much 
later recognized that the conditions for MI would be significantly modified for discrete settings relevant to, for instance, 
the local denaturation of DNA [5] or coupled arrays of optical waveguides [6] . In the latter case, the relevant model 
is the discrete nonlinear Schrodinger equation (DNLS), and its MI conditions were discussed in [7]. Most recently, 
the MI has been recognized as responsible for dephasing and localization phenomena in the context of Bose-Einstein 
condensates (BEC) in the presence of an optical lattice [8-11]. 

In this brief report, we present an alternative approach to the modulational stability of plane waves in the context 
of the nonlinear Schrodinger equation 



g ■ ii> t = -V** - u\m, (i) 

^-j- , where ip is a complex field, the subscripts denote partial derivatives with respect to the corresponding variable and 
• U is a constant prefactor (the strength of the nonlinearity) . We examine MI using a time dependent variational 
approach (TDVA), and study the results in comparison with the standard linear stability (LS) calculations. It should 
be mentioned that the use of TDVA for the study of solitons at the classical [12] and even at the quantum [13] 
level is not novel. What distinguishes our study from these earlier ones is the use of the Mi-motivated ansatz in the 
TDVA (see below). We also note in passing that Ml and solitons from a quantum mechanical point of view have 
been considered in a number of references including (but not limited to) [14]. Furthermore, a similar in spirit, 3-mode 
approximation was systematically developed in the works of [15]. There are however a number of differences between 
the latter and the present approach such as e.g., our use of the variational formulation of the problem (instead of the 
application of the 3-mode ansatz in the dynamical equation in the context of [15]), as well as the fact that we are 
perturbing around an exact plane wave solution, while in the case of [15], the plane wave is an additional mode in the 
relevant expansion. 

In the linear stability framework (see e.g., [7,9] and references therein for relevant details), the stability of the plane 
' waves has been examined. The latter are of the form 

tp(x, t) — ipo exp [i (kx — cot)] , (2) 

and constitute exact solutions of the nonlinear Schrodinger equation with a dispersion relation 

lo = k 2 - Uipl. (3) 

Then, the MI is examined in the LS framework using the linearization 

u(x, t) = (ipo + ec) exp[i((kx — ut) + ed(x, t))] (4) 

and analyzing the 0(e) terms as 
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c(x,t) — c exp(z/3(x, t)), d(x,t) = d cxp(i[3(x,t)). (5) 

Using f3(x, t) = qx — fit, the dispersion relation connecting the wavenumber q and frequency fl of the perturbation 
(see e.g., [1]) 

(-fl + 2kq) 2 = q 2 (q 2 - 2Ulft) (6) 

is obtained. This implies that the instability region for Eq. (1) appears for perturbation wavenumbers q 2 < 2Utp 2 , 
and in particular only for focusing nonlinearities (to which we restrict this study). 

We now attempt to identify the interval of unstable wavenumbers by means of the TDVA. In particular, we start 
from the Lagrangian L 



L 



i ( ^-M*)-iVxi 2 + |ivi 4 



dx, (7) 



and consider a modulation of the plane wave of the form 

i> — H>o + a (t) ex p(*0o(i)) exp(iqx) + b(t) exp(i(^,(i)) exp(— iqx)] exp (i(kx — cot)) . (8) 

However, instead of considering the modulation directly at the level of the equation, the variation of our approach is 
that we use the modulational ansatz in the Lagrangian. This constitutes the basic novel ingredient of this variational- 
type approach to the modulational instability. 

Here we consider an annular (1-dimensional) geometry, which imposes periodic boundary conditions on the wave- 
function ip(x) and integration limits < x < 2tt in Eq. (7). This results in the quantization of the wavenumbers 
k,q = 0, ±1,±2, .... However, it is clear that our results can be easily generalized to the case of an infinite, open 
system and to higher dimensions. 

After substitution of Eq. (8) into Eq. (7), we obtain the variational Lagrangian 

L = Tr[-2{a 2 ^ a + b 2 <j> b ) + 2{UiPl - q 2 ){a 2 + b 2 ) - - 4qk(a 2 - b 2 ) + 4U^abcos(4> a + fa) 

+U(a 4 + b 4 + 4a 2 b 2 )]. (9) 

It is clear from this Lagrangian that the pair <p a (t) , 4>b{t) can be interpreted as the generalized coordinates of the 
system, while A(t) = 2a 2 (t), B(t) — 2b 2 (t) are the corresponding momenta. In particular, the pairs A(t) , <j> a {t) and 
B(t),(f>b(t) are canonically conjugate with respect to the effective Hamiltonian 

H eff = (UiP 2 - q 2 )(A + B) — 2qk(A — B) + 2U^VABcos((f> a + <j> h ) + j(A 2 + B 2 + AAB), (10) 

which is an exact integral of motion on the subspace spanned by Eq. (8). 
The Lagrangian equations of motion are: 

a<j) a = C\a + C 2 bcos(<p a + <f> b ) + Ua(a 2 + 2b 2 ) (11) 

>a = C 2 bsm{<j) a + (j)b) (12) 

b^ b = C 3 b + C 2 a cos(^ a + 4> b ) + Ub(b 2 + 2a 2 ) (13) 

> b = C 2 asm((f) a + <j> b ), (14) 

where C\ — Uip 2 — q 2 — 2qk, C 2 = Uip 2 , and C3 = Uip 2 — q 2 + 2qk are constant prefactors. 

If we now keep all terms to 0(a) in Eqs. (11)-(14) [which is consistent with an approximation linear in a], we obtain 

a = b (15) 
a = C 2 asm((j)) (16) 
0=(C 1 +C 3 ) + 2C 2 cos(0), (17) 

where cf) — 4> a + 4>b- The latter equation has the solution 
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4>(t) = 2 arctan 



V Wg - g 2 ) g 2 



tanh U(2Uxp 2 - q 2 )q 2 t 



(18) 



Two different cases arise here, corresponding respectively to whether the instability criterion is satisfied or not. 
Namely, when 2UipQ — q 2 < the solution of Eq. (16) is 



i(t) ~ sjl Sin (7 (^ - 2UiPl)q*tj, 



(19) 



while when 2Uipo — Q 2 > the solution of Eq.(16) is 



a(t) ~ Jl + ^ sinh ( v /(2[/^o 2 - <7 2 )<Z 



(20) 

The solutions signal the appearance of the modulational instability when the threshold condition q 2 = 2UipQ is 
crossed (passing from higher to lower perturbation wavenumbers) . This is also clearly shown in the time evolution of 

a(t) in accordance with Eqs. (19)-(20) also shown in Fig. 1 in the case of q = 2 for 2U q f = 0.2 (see the left panel of 
Fig. 1) and 2U q f° = 1.2 (see the right panel of Fig. 1). 




FIG. 1. The left panel shows the (stable oscillatory) time evolution of a(t) for the case 2U f a = 0.2, in accordance with 



Eq. (19). The right panel shows the (unstable, exponentially growing) time evolution of a(t) for the case of 2U f° = 1.2, in 
accordance with Eq. (20). 



An alternative, more intuitive way to appreciate the linear stability result from a dynamical systems viewpoint. 
This consists of reducing Eqs. (11)-(14) to a one degree of freedom setting with an effective potential energy landscape 
whose (parametric) variation will elucidate the instability. Along these lines, using A(t = 0) = in (without loss of 
generality) in Eq. (9) and A(t) = B{t) (from Eqs. (12) and (14)) in Eq. (10), we have: 



H = 7T 



2(f7Vo - 1 2 ) A + 2 UA2 + 2^0^ cos < 



= 0. 



(21) 



Eliminitating <j> from Eqs. (12) and (21) for A(0) = 0, we obtain the "energy equation" for A 

\A 2 + V ef f = 0, 



(22) 



where the effective potential V e f / (A) is of the form: 



V eff (A) = 2<z V - 2U^)A 2 + 3C/(C/Vo - Q 2 )A 6 + ^U 2 A 4 



(23) 



One can then examine the stability of the effective potential by evaluating its curvature at A = 0. We thus obtain 
V"ff(A)\A=o = 4<7 2 (<j 2 — 2UiPq) and hence the potential will be convex (and therefore the dynamics will be stable) 
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for q 2 > 2UipQ, while it will be concave (and the dynamics unstable) for q 2 < 2Utp 2 . Hence in this case also, we 
retrieve the modulational stability criterion. The effective potential is shown for the modulationally stable, unstable 
and marginal case in Fig. 2. 




FIG. 2. The effective potential of Eq. (23) is shown for as a function of A for U = tpo = 1 and three different values of 
q: q — 1 (modulationally unstable; solid line), q — \/2 (at the threshold; dashed line) and q = 2 (modulationally stable; 
dash-dotted line). 

Now we turn to a more interesting case, where the coefficient of the dispersion term, as well as the coefficient of 
the nonlinear term in Eq. (1) are temporally modulated, namely we examine the equation 



i^Pt = -D{t)1> xx - u(tM 2 ip. 



(24) 



Our aim is to derive the modulational stability equation via the TDVA, for general D(t) and U(t). It is interesting 
to note that this equation has become of increasing importance in the past decade due to applications both in 
optics and also, more recently, in soft condensed-matter physics. In particular, in optics, the case of D(t) periodic 
and U(t) constant is of relevance in the context of the so-called dispersion management. The latter is based on 
periodic alternation of fibers with opposite signs of the group- velocity dispersion [16]. We note in passing that in this 
application time t is, in reality, the propagation distance (i.e., space), while x corresponds to a retarded time variable. 
An alternative setting where D(t) is constant, but U(t) can be temporally modulated (via a Fcshbach resonance, i.e., 
an external magnetic field; see e.g., [17]) can be found in Bose-Einstein condensation. In the latter setting, there has 
been an explosion of interest recently in time dependent scattering length and its effect on patterns, coherent structures 
and collapse thereof (a number of very recent references can be found in [18-21]. While our primary motivation in 
considering MI through the TDVA in Eq. (24) principally stems from this recently explored experimental potential 
in Bose-Einstein condensates, we should note that this type of problem was investigated earlier in nonlinear optics, 
see e.g., [22]. 

We consider the perturbation of the form 



if) = ipi [1 + w(t) cos(gx)] 



(25) 



to the plane wave solution ipi = k fo D ( s ) ds +/ u(s)ds+kx) ^ ^ Q ^ ce ^hat in Eq. (25), we are using for simplicity a 
variant of the ansatz of Eq. (8), with a — b, 4>a — 4>b an d (without loss of generality) ipo = 1- Following the same 
procedure as above, we can obtain the stability equations for w = w r + iwf. 



w r = q 2 D{t) Wl ~-U{t) (wf 



- Wiw; 



Wi = - [q 2 D(t) - 2U(t)] w r + -U(t) (w^ + w r w 2 ) 
Hence, at the linear level, we can derive the following stability equation 



Pit) 
D(t) 



w r - q 2 D{t) [q 2 D(t) - 2U(t)] w r 



(26) 
(27) 

(28) 



By determining the windows of stability of the ordinary differential equation of (28) , the modulational stability of Eq. 
(24) is determined. It is further worth noting that for D{t) constant and U{t) time-periodic, Eq. (28) falls into Eq. 
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(2) of [21] and becomes Hill's equation for which many stability results are known in the mathematical literature [23]. 
Furthermore, in the case of U(t) = 1 + 2acos(wi), Eq. (28) falls into Eq. (2) of [18] and the resulting equation is of 
the Mathieu type for which explicit stability windows can be computed (for details see [18] and references therein). 
It then becomes naturally an interesting problem in mathematical physics to determine the stability of Eq. (28) for 
more general cases (e.g., with both coefficients periodically varying etc.). 

In this brief report, we have revisited the modulational instability from a different point of view, namely a vari- 
ational one. We have used this dynamical systems' type approach to derive the Euler-Lagrange equations for the 
time-dependent perturbation ansatz parameters and have examined their stability for different wavenumbers of the 
perturbation (which affect the constants of the ensuing set of ordinary differential equations). We have retrieved, 
in a simple and intuitive way, the criterion for the instability. The technique has also been generalized in cases in 
which the coefficients of the dispersion and/or nonlinearity are temporally varying (a case which we have argued to be 
relevant to a variety of applications). We have found the corresponding stability condition obtaining a novel ordinary 
differential equation, whose special cases correspond to stability/instability criteria established previously It would 
be interesting to extend the considerations of this method (which seems applicable to any setting with an underlying 
Lagrangian/Hamiltonian structure) to contexts with explicit spatial dependence of the potential (see e.g., [9,11]). 

The support of NSF (DMS-0204585), UMass and the Clay Institute (PGK) is gratefully acknowledged. Work at 
Los Alamos is supported by the US DoE. 
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